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Abstract 

Double scaling may be observed in the Drell-Yan process at the Tevatron and in the 
AblCE detector at the Large Hadron Collider in a window of masses. In the double 
scaling limit of the cross section, the higher order QCD corrections are quite large, and 
are driven by the rise in gluon densities in the double scaling regime. The naive parton 
model cannot be valid even at asymptotic energies in the region of phase space where 
double scaling holds. 
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The HERA experiments have observed double asymptotic scaling M of the deep inelastic 
structure function F 2 (x, Q 2 ) by probing a region of kinematics with large momentum transfer, 
Q 2 , and small Bjorken variable, x. In this paper, we discuss the possibility of future hadron 
colliders observing similar double scaling phenomena. We examine the simplest such process 
here — that of Drell-Yan production of a lepton pair. The appropriate kinematic range will be 
accessible to the ALICE experiment at the CERN Large Hadron Collider and to experiments 
at the Tevatron. 

In this process a lepton pair is produced in a hadron collider at a centre of mass energy 
\/S. The kinematics of the pair is specified by its invariant mass, Q, and the longitudinal 
momentum fraction, x F = 2q z /\/~S. If x F = 0, then at large \/S and Q, the ratio r = Q 2 / S 
can be small. The process will probe parton densities in the double scaling region x\ = x 2 = 
y/r <C 1 and scale Q 2 large. The pair mass must satisfy the conditions Qo < Q < xoVS 
with appropriate choices of Qq and Xq. In this region the higher order QCD corrections are 
computable in the double asymptotic scaling (DAS) ]I[ approximation and are not negligible. 
This is what we show in this communication. 

The Drell-Yan cross section is the sum of the qq annihilation process and the QCD 
Compton process ||. The cross section at leading order (LO) in QCD [i.e up to 0(a s ) in 
the strong coupling constant) comes from the diagrams shown in Figure [TJ. We treat the two 
sets separately. Take the double differential cross section for the annihilation contribution 
to the Drell-Yan cross section — 
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The DAS approximation must be consistently applied to the whole expression, i.e., to both 
the "finite parts" , 
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and the parton density combination, 



H(xi, x 2 , Q 2 ) = e / [lf( x i^ Q 2 )Qf( x 2, Q 2 ) + (xi <-> x 2 ) 
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Before moving on to technicalities, we would like to set out the justification for this claim. 
Recall that in defining the parton densities through the DIS process, we absorbed collinear 
divergent parts of the higher order correction into the definition of these densities through 
renormalisation group equations which are the DGLAP equations. Of the collinear singular 
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part, the terms proportional to log(<5 2 /A 2 ) give the splitting functions (Mellin transforms of 
the anomalous dimensions of the leading twist operators) and the remainder are called "finite 
parts". These latter are process dependent. In the DAS approximation for deep-inelastic 
lepton scattering, the DGLAP equations are solved analytically by expanding the anomalous 
dimensions about their right-most singularity at n = 1 in Mellin space i.e they are expanded 
in Mellin space in the quantity A = n — 1. The leading terms are of the form 1/A (and come 
entirely from the splitting functions), and the first sub- leading terms are independent of A. 
In the Drell-Yan cross section in eq. (||) some of the collinear divergent parts of the process 
have been absorbed into the scale dependent parton density combination in eq. (^) . If this is 
treated in the DAS limit, then certainly the remaining part of the collinear divergent terms, 
in eq. (§), must be given the same treatment, to wit the 1/A and constant terms must be 
retained in the Mellin transform (with respect to r) of eq. ([!]). 

In order to be able to use the DAS forms of the quark and gluon distribution functions 
in the double probability distributions given in eq. (^), we need to first express the quark 
distributions in terms of the singlet (S) and non-singlet (NS) components. We define these 
through the equations 
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Here n F is the number of active flavours. This allows us to express the quark density for 
flavour /, qj, in terms of q s and qnsj- In the DAS limit where the non-singlet can be 
neglected (up to sub leading order) compared to the singlet — qj is just given by the singlet 
quark densities. 

Thus in the DAS limit the DGLAP equations can be solved exactly to give 



Q(x,Q 2 ) = xg(x,Q 2 ) - 
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Q(x,Q 2 ) = xq s (x,Q 2 ) = N'g(x,Q 2 h/p. (7) 

where pa = log(xo/x) and a/p = log(t/to) with t = ln(Q 2 /A 2 ) (t which is t evaluated at 
Qq, and Xo are chosen so that the approximation holds). The usual gluon density function 
is denoted by g(x, Q 2 ). The constants appearing in these equations are given by 
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and the running coupling a s (t) = Air / (flat) . N and N' are overall normalisation constants 
which will be fitted with data or parametrisations. These are expressions valid to lowest 
order in the strong coupling constant g s . 

There are no 1/A terms in the Mellin transform of f q and f q . The constant terms are 
identifiable by taking the first moment of these quantities — 

M[f q ] = - 7 -+0(A), and M[f q ] = -1 + 0(A). (9) 
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As a result, the leading part of the finite terms can be modelled by the z-space functional 
forms — 

f q ( z )~J-5(l-z) and f q (z,y*)~-6{l-z). (10) 

The delta-functions tell us that this subleading correction is due entirely to the soft part of 
the phase space in the QCD corrections. 

Substituting these into eq. ([]]) then gives us — 
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The appropriate expression for the parton density combination is 
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where x\ = x 2 = y/r, and hence cri = o 2 = \og(x / y/r). The overall normalisation constant 
Hq can be fixed from data on parton densities. 

We turn next to the contribution of the Compton process to the cross section. The full 
expression is — 
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and the parton density combination, 

G(x 1: x 2 , Q 2 ) = g(x 2 , Q 2 ) e / [9/(^1, Q 2 ) + Qf( x u Q 2 ) 

f 

are both to be evaluated consistently in the DAS limit. 

The leading parts of the finite terms are identifiable from the first moment- 
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M[f g } = -- + 0(A), 



and M [f 9 } = ~(l + y* + ^j+0(A). 
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In z-space, they can be modelled up to subleading order in A by the functions 
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Substituting into eq. ( p"3D gives the DAS result- 
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The appropriate expression for the parton density combination is 
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The constant Go can again be fixed from the numerical values of the parton densities. 
This concludes the computation, since the Drell-Yan cross section 
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and the two terms on the right are given in eqs. (11,19). We make three remarks- 
• Since x F = 0, we have x\ = x\ = y/r, and hence 
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o"i = o"2 = vlog(t/to) log^o/v 7 ^)- 
Also, t = \og(Sr / A 2 ) . Hence, the only free variable in eqs. ([11,191) is r. 
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• In general, QCD corrections involve parton densities at all momentum fractions in the 
range y/r < x < 1, whereas in the DAS only the densities at x — yfr are involved. This 
is due to the delta-functions in eqs. (pTD| , |T8D . The same physics reduces the DGLAP 
equations for structure functions to the wave equation (T|]. 

• The coefficient of a s in eq. (|TT|) is negative and larger than unity. This is not a problem, 
since this cross section is not physical, being only one part of the total physical cross 
section. Adding eq. (|T9|) makes the total physical cross section positive. 

Since we examine DAS at leading loop order, the parton densities to be used in fixing the 
normalisations Hq and Go must be obtained at LO. For this reason we work with the GRV 
94 LO pi set of parton density parametrisations. We find that DAS is valid if we choose 
Ql = 1 GeV 2 and x = 0.1. It turns out that G /H ps 40. Such a large value of this ratio 
is generic, since the gluon densities are substantially larger than the quark densities. As a 
result, the Compton process dominates over the annihilation process in the DAS. 

Recall that the increase of the Drell-Yan K-factor over unity is due to the parts of the 
cross section proportional to a s . From the DAS expressions, and the fact that G 3> H , it is 
clear that the K-factor is much larger than unity for — > oo. In fact, from the behaviour 
of a s , one can see that the K-factor vanishes as \ jt. The rest of the factors in eqs. ( |n| , |i~9"D are 
almost independent of S at fixed r. Since the K-factor in the DAS remains non-negligible at 
the largest accessible energies, the naive parton model cannot be recovered wherever DAS 
is valid. In contrast, the K-factor at larger values of r are rather close to unity at these 
energies M, and the parton model may be taken as a reasonable approximation to reality. 
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In Figure |2| we compare the DAS forms of the cross section with the full one-loop expres- 
sion. It is interesting to note the disagreement between the two sets of computations at low 
r. This is expected, since low r corresponds to low Q 2 , where DAS breaks down. Similarly, 
at large r, DAS breaks down because x is not small enough. Thus, the applicability of DAS 
to Drell-Yan is bounded both above and below, in r. This is to be contrasted with the 
situation in DIS. 

In Figure [5] we show the differential cross section Sda / dlogrdxp as a function of y/r for 
three different values of \/S. Note that the cross section for these three values of \/S lie 
almost on one single universal curve. This is the double scaling curve. Departures from DAS 
are visible at \/t ~ x and rS ~ Q^. 

In the DAS limit of hadronic cross sections, the leading correction is universal — this is 
the 1/A piece which comes entirely from the splitting functions, i. e., the collinear region 
of phase space, and is absorbed into the parton densities. The next term in the expansion 
in A is not universal. However, since it arises in the soft part of phase space, it should be 
easier to evaluate and sum the DAS phase space to all loop orders than to evaluate the full 
correction. In fact, this resummation might even be necessary. The reason is in the large 
value of the K-factor, as mentioned above. Since the one loop correction is so large, it may 
actually be necessary to sum the logQ 2 log(l/r) contributions to all orders, so as to describe 
the physical process. This is an extension to hadronic collision of the usual DLL procedure 
in lepton-proton DIS. We leave this for the future. 

One of us (RB) would like to thank D. Indumathi for useful discussions. This work was 
started at the 5th Workshop on High Energy Physics Phenomenology (WHEPP-5), held in 
January 1998 at Pune, India. 
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Figure 1: The diagrams that contribute to the Drell-Yan process up to 0(a s ) in the strong 
coupling constant — (a) the annihilation diagrams and (b) the Compton diagrams. 




Figure 2: The DAS prediction for the Drell-Yan cross section as a function of y/r for x F = 
compared with the full cross section at three different values of \/~S. In each case, the DAS 
curve is the lower one. The curves for 16 TeV and 1.8 TeV have been scaled by a factor of 
100 and .01 respectively for clarity in display. 
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Figure 3: The DAS prediction for the Drell-Yan cross section as a function of yfr for x F = 
at the values of v^S* shown. 
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